Nonlocal spin correlations mediated by a superconductor 
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Nonlocal charge correlations induced in two normal metals contacted separately to a supercon- 
ductor have been studied intensively in the past few years. Here we investigate nonlocal correlations 
induced by the transfer of pure spin currents through a superconductor on a scale comparable to 
the superconducting coherence length. As with charge currents, two processes contribute to the 
nonlocal spin signal: crossed Andreev reflection (CAR), where an electron with spin up injected 
from one normal metal into the superconductor results in a hole with spin down being injected into 
the second normal metal, and elastic cotunneling (EC), where the electron with spin up injected 
from the first normal metal results in an electron with spin up being injected into the second normal 
metal. Unlike charge currents, however, the spin currents associated with CAR and EC add due to 
the fact that the bulk superconductor cannot sustain a net spin current. 

PACS numbers: 74.45. +c 72.25.Ba 72.25.Hg 72.25. Mk 



Electrons in two spatially separated normal metals in 
contact with a superconductor show nonlocal correlations 
that are mediated by their mutual interaction with the 
superconductor [flS]. Two processes are responsible for 
these correlations. In the first process, called crossed An- 
dreev reflection (CAR), an electron with one spin orienta- 
tion, e.g., spin up, is injected from the first normal metal 
(N) into the superconductor (S), resulting in a hole with 
spin down being injected from the superconductor into 
the second normal metal, with a concomitant generation 
of a Cooper pair in the superconductor. In the second 
process, elastic cotunneling (EC), the spin- up electron 
injected from the first normal metal gives rise to a spin- 
up electron being injected into the second normal metal. 
Both processes are exponentially suppressed with the dis- 
tance between the two NS interfaces on a length scale of 
the order of the superconducting coherence length £g [5] . 
In the case of charge currents, the net current injected 
into the second normal metal in response to the drive cur- 
rent injected from the first normal metal is the difference 
of the two contributions CAR and EC. As the relative 
amplitude of CAR and EC is predicted to depend on the 
transparency of the interface [5J, the effect of electron- 
electron interactions [7] , and so forth, the sign of the net 
current injected into the second normal metal may be 
positive or negative. Experimentally, nonlocal correla- 
tions due to CAR and EC in charge transport have been 
verified by many different groups [8HTU]. 

In this paper, we study theoretically nonlocal correla- 
tions induced in a NSN structure in response to a pure 
spin current. As with charge currents, both CAR and EC 
contribute to the resulting nonlocal spin signal. However, 
due to the fact that one cannot inject a pure spin current 
into the bulk of a s-wave superconductor (which is the 
case of interest here), we find that the resulting nonlocal 
spin current is the sum of CAR and EC contributions: 
Injecting a spin-down hole (CAR) or injecting a spin-up 
electron (EC) into the second normal metal, both corre- 
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FIG. 1: Schematic diagram of the device. The bias current 
is applied from a ferromagnet Fi to a normal metal Ni and 
nonlocal voltage is measured between F2 and N2. 



spond to injecting a net up spin. Consequently, with a 
combination of charge and spin transport measurements 
on the same device, one should in principle be able to 
separate the contributions of CAR and EC, which has 
not been possible to do so far. 

In addition to CAR and EC, there are additional pro- 
cesses that may contribute to the nonlocal signals in NSN 
or FSF structures. These are charge [TTHlT)] and spin 
[mill] imbalance, associated with the injection of quasi- 
particles, with energies larger than the superconducting 
gap , into the superconductor. In contrast to these stud- 
ies, the proposed experiment addresses the regime of sub- 
gap transport. Thus, no quasiparticle is injected above 
the gap, and the effect does not depend on the long spin 
relaxation times recently observed in the nonlocal mea- 
surements of Refs. [HI [17] at applied bias voltage larger 
than A/e. 

Figure IT] shows the schematic view of the device geom- 
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FIG. 2: Nonlocal correlations near N1SN2 structure. Pure 
spin current is injected from Ni and a finite spin current is 
generated in N2 by CAR and EC. One should note that the 
direction of the spin currents attributed to CAR and EC are 
in the same direction and there is no net spin accumulation 
in the superconductor. 



etry that we consider here. Our device is similar to the 
usual NSN hybrid structure, except that there are two ad- 
ditional ferromagnets, Fi and F2, placed in contact with 
the two normal metals, Ni and N2, respectively. One of 
the normal-metal/ferromagnet junctions (N1F1) serves 
as a spin injector. If one drives a bias charge current / 
from Fi to Ni, as shown in Fig. [TJ spin accumulates at 
the N1F1 interface. The resulting spin imbalance leads 
to spin diffusion in all possible directions in Ni. Since 
the charge current is drained from the left side of Ni, 
the right side of Nj carries a pure spin current with no 
net flow of charge. The second NF interface (N2F2) is 
used to detect the spin current flowing through N2 by 
measuring the nonlocal voltage V n \ that establishes it- 
self between N2 and F2. If Fi and F2 are designed to 
have different coercive fields, it is possible to realize both 
parallel and antiparallel magnetization directions by ap- 
plying an external magnetic field. Below we determine 
the nonlocal spin signal by computing the nonlocal resis- 
tance R n \ — V n i/I. 

In the following, we assume that the normal metals are 
oriented along the cc-axis whereas the ferromagnets are 
oriented along the y-axis (see Figure pi). To simplify the 
notation, we use two different coordinate systems for Ni 
and N2 with origins at the respective NF interface and 
the x-axis directed toward the superconductor, i.e., the 
respective NS interfaces are at x = Lj (i = 1, 2). 

Let us first consider the processes that occur at the 
N1SN2 interfaces in more detail. As shown in Fig. [2] 
an incoming electron with subgap energy may undergo 
four possible processes at the NiS interface: local An- 
dreev reflection (AR), normal reflection (NR), crossed 
Andreev reflection (CAR), and elastic cotunneling (EC). 
Only the nonlocal processes (CAR and EC) contribute 



to the nonlocal signal. If the length of the supercon- 
ducting link between Ni and N2 is much longer than £g, 
these processes are suppressed. In that limit, only AR 
and NR occur at the NiS interface, and thus the detec- 
tor voltage between N 2 and F 2 is zero. In the opposite 
regime, when d is shorter than £<j, however, nonlocal pro- 
cesses such as CAR and EC can occur and generate a 
spin current in N2. The resulting spin accumulation at 
the N 2 F 2 interface yields a finite voltage difference be- 
tween N 2 and F 2 . It is important to note here that, in 
the subgap regime, the superconductor cannot support 
spin accumulation: the spin of an electron injected into 
the superconductor at the NiS interface is transferred to 
N2, either via a hole with the opposite spin (CAR) or via 
an electron with the same spin (EC) , thus leaving no net 
spin in the superconductor. The spin current injected 
from Ni into S, Ifc s , equals the spin current injected 

from S into N 2 , — ^n 2 S- Here J§. s = iJ-S — is * s ^ ne 
difference of spin-up and spin-down currents, with a fixed 
spin quantization axis collinear with the magnetizations 
in F]/ 2 . As a consequence, as can be seen schematically 
in Fig. [2j both CAR and EC contribute to a spin current 
in N2 in the same direction, such that their contribu- 
tions to the spin signal add. This is completely different 
from the case of charge current, where one measures the 
difference between CAR and EC. Overall, the spin cur- 
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rent through the superconductor, 

is related to the difference between the spin imbalances 

S/iNiiLi) = [//Jj.(.Lj) — /Ujj.(Lj)]/2 at the two interfaces, 



Ii = ^-[S^ALi) 



6h-n 2 (L 2 )], 
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Here /% are the spin-resolved electrochemical poten- 
tials. Note that /XN 1 (ii) = fJ,-^ 2 (L 2 ), with /j,^. (Li) = 
[/zjj. (Li) +/%. (Lj)]/2, in the absence of a charge current. 

Furthermore, Gg = Gcar + Gec, where Gcar/ec are 
the conductances due to CAR and EC, respectively. 

To describe spin accumulation and spin transport at 
the NF interfaces, we use the model introduced by Taka- 
hashi et al. 18 . The current density j° = — (cr° / e)\7 fi a 
for spin a = j" , 4, is related with the spin resolved conduc- 
tivity a^, and electrochemical potential \i a a in material 
a = Fj,Nj. The continuity equation for the charge cur- 
rent density j a = j a + j a reads 



Vj Q = 0, 



(2) 



while the relaxation of spin imbalance is described by the 
phenomenological equation 



V 2 <W - (l/X a )5^ a = 0, 



(3) 



where X a is the spin diffusion length. The spin current 
density is j a = j£- j a . 



The general solution of Eqs. (pi) and (pf in Fj reads 
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(4) 



with I\ = I and J2 = 0. 0(y) is the Heaviside function. 

Note that the chemical potentials in the ferromagnets 

are defined in the same way as in the normal metals, i.e., 

1*/ 
Mf (v) — ^Fi{y) i SfiFiiy)- For simplicity, we assume 

identical cross-sections Af, spin-diffusion lengths Af, and 

i* i 

total conductivities op = ct f . + er^. in the wires Fi and 

F2, while pp. = (crp. — ct f .)/ctf. In particular, p Fl and 
Pf 2 have the same modulus pp , while their sign is deter- 
mined by the magnetization direction in Fj, p Fi = ±pp, 
depending on whether electrons with spin up/down cor- 
respond to majority/minority (+) or minority/majority 
(— ) electrons, respectively. Using Eq. Q to compute 
the current densities at y = O*, we can relate the spin 
imbalance at the FN interfaces to the drive current and 
the spin currents I F . N . injected from Fj to N* through 
the F,N,- interface, 



<5^ F ,(0) = 2 e - Rp (PFji ~ ^nJ 



(5) 



where Rf = Af/[Afof(1 — P f )]- 

The currents through the F^N^ interface may also be 
expressed in terms of the potential drops across the in- 
terface, 
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(MFiNU + PT i S/J,F i Ni) > (6a) 

(PTi^FiN, + SfJ.FiNt) , (6b) 



where G T , = G^+G^ andp Ti = (G^-G^)/G Ti are re- 

T / I 

lated with the (tunnel) conductances G^ of the N^F^ in- 
terface for spin up and spin down electrons. Furthermore, 
MFiNi = Mf, (0) - ft, (0) and(5^ FjNi = Sfi Fi (0) -<fyt Ni (0). 
The sign of pr t again depends on the magnetization di- 
rection of the ferromagnet. Note that spin- flip scattering 
at the interface has been neglected in Eqs. (pi) for sim- 
plicity. Inverting these equations, we obtain 

^FiN, = eR Ti (Ii ~ PTjFiNi) > ( 7a ) 

^FiN, = eR Tl {-Prji + ^f,nJ > ( 7b ) 

where i? T , = l/[G Ti (l -p|.)]. 

The general solution of Eqs. (pi) and (pi) in the normal 
metals Nj reads 

g J T 

MN 4 (a;) = MNi(O) -j— 9(-x), (8a) 

1 / _M l^i-^l \ 

^ Ni (or) = -ei? N M Fl N s e ^-^N^e An J ( 8b ) 



with i?N = 1/Gn = An/cn^n- As for the ferromagnets, 
here we also assume identical cross-sections A^ and spin- 
diffusion lengths An, and conductivities ctn in the wires 
Ni and N 2 . Furthermore, o" N = cjj = ctn/2- Using these 
equations, we find in particular 

SmM = \eR N {l ¥ini -I^ Si e- Li '^), (9a) 

Sm(Li) = ^ei? N (/# 4Nj c- £ «/ An - 7 NiSi ) . (9b) 

Combining Eqs. (p| with (IT]), we can eliminate the spin 
currents and imbalances at the NS interfaces to obtain 
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2(G++G N ) 

„-(L 1 +L 2 )/\ l . 



(10) 



where we used the notations 1 = 2 and 2 = 1. 

The nonlocal voltage V n i between F2 and N2 is given 
by /jn 2 (+oo) — Mf 2 (— 00) • Using Eqs. Q and (J8J), we find 
Uni = [mf 2 n 2 +PF 2 ^/ i F 2 (0)]/e which, using Eqs. (pi) and 
( pa) , may be rewritten as V n i = -(pF 2 /pF)-R NFi / F2N2 
where we defined i?NF; = Pf-Rf/2 + \pT t \Rr t and used 
PtJ\pt 1 \ = Pf s /pf- Finally, combining Eqs. J5}, ( pb| , 



and (10), we determine the spin current I F N as a func- 



tion of the injection current / to obtain the nonlocal spin 
resistance, 
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±- 



i? NFl fl N F 2 -Roe- (Ll+£2)/AN 

^NSi-RnS 2 



^2 e -2(L 1 +L 2 )/A N 



(11) 



where R = G+/[4G N (G+ + G N )] and R NSi = i? N /2 + 
Rf /2+R T% -R er 2L '/ XN . The overall sign of R n \ depends 
on whether the ferromagnets are aligned parallel (+) or 
antiparallel (— ). 



Equation ( 11 ) is the central result of this paper: it pre- 



dicts a finite nonlocal resistance with no charge current 
injected into the superconductor. The exponential de- 
pendence ex e _ ( Ll+L2 )/ AN of i? n i on the lengths of the nor- 
mal wires clearly shows that it is due to the spin transport 
through the structure. In the case of fully polarized fer- 
romagnets, where \pT t \ = Pf = 1 and -Rf/2 + R^ -^ 00, 
the result (pTj) simplifies to R n \ = ±i? e~ (Ll+L2)/AN . 



A nonloncal spin signal is possible also in the absence 
of superconductivity. Note that, at Gg 3> Gn, when the 
contribution of the superconducting element to the non- 



local signal is negligible, Eq. ( 11 1 is similar to the result of 
Ref. [TJ5] up to factors due to a different geometry of the 
normal part of the FNF spin valve. In our geometry, a 
conventional spin valve is realized at temperatures above 
the superconducting transition temperature T c . We ex- 
pect the decay lengths of the nonlocal signal within the 
superconductor to be quite different below and above the 
transition. At T <C T c , our results yields R n \ oc e~ d ^ s 
since Gcar/ec decay exponentially on the scale £g. In 



contrast, above T c , the nonlocal resistance should be pro- 
portional to e~ d ' Xs , where As is the spin diffusion length 
of the superconductor in the normal state. Typically, As 
for Al is - 500 - 1000 nm and f s is ~ 100 - 300 nm. 
Thus, one would expect an abrupt change in the nonlo- 
cal resistance when the superconductor transitions to its 
zero-resistance state. This effect is possible to demon- 
strate in experiments similar to Ref. [8]. 

The magnitude of the nonlocal resistance that is in- 
duced by the injection of pure spin current may be es- 
timated from the nonlocal signal measured in the case 
of charge current injection. In that case, based on the 
formalism suggested by Falci et al. [3], the nonlocal re- 
sistance can be written as i?^ SN = (Gec — Gcar)/G ar , 
where the conductance due to AR at a single NS inter- 
face Gar ^> Gcar, Gec is assumed to be the same for 
the NiS and N 2 S interfaces. Using the measured i?J^ SN as 
well as estimated values of Gar from Cadden-Zimansky 
et al. [TO] yields a rough estimate of Gec — Gcar ~ 
0.5 f2 _1 . A lower bound for the nonlocal spin resistance 
can be obtained if we assume that Gec ^ Gcar- For 
copper wires with a spin diffusion length An ~ 1 £tm 
and cross section A-^ =100x50 nm 2 , corresponding to 
Gn ~ 0.3 fi -1 , the factor Rq in Eq. (11 1 would be of the 
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order of 0.5 il in this case. If Gcar and Gec are of the 
same order of magnitude, the nonlocal resistance would 
likely be much larger. 

In summary, we have investigated the nonlocal signal 
that arises from the injection of pure spin current into a 
superconductor. We have shown that a finite electrical 
resistance may be generated entirely due to the nonlocal 
correlations mediated by a superconductor. In contrast 
to other recent work, this nonlocal signal arises from spin 
transport at energies far below the superconducting gap. 
Measurements of the nonlocal resistance resulting from 
charge and spin injection on the same sample would allow 
independent determination of the contributions due to 
CAR and EC. 
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